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Abstract 

Using a kinematic approach we show that the non-adiabatic, non- 
cychc, geometric phase corresponding to the radiation emitted by a three 
level cascade system provides a sensitive diagnostic tool for determining 
the entanglement properties of the two modes of radiation. The nonuni- 
tary, noncyclic path in the state space may be realized through the same 
control parameters which control the purity /mixedness and entanglement. 
We show analytically that the geometric phase is related to concurrence 
in certain region of the parameter space. We further show that the rate of 
change of the geometric phase reveals its resilience to fluctuations only for 
pure Bell type states. Lastly, the derivative of the geometric phase carries 
information on both purity /mixedness and entanglement /separability. 

1 Introduction 

Geometric phase [T] and topological phases have provided impetus[2] to 
quantum information processing. Of the two geometric phase is easier for 
implementation. Both these approaches may be designed to realize fault 
tolerance in addition to resilience to decoherence [H [5] . Simple quan- 
tum gates using geometric phase have been demonstrated experimentally 
in the nuclear magnetic resonance set up [6]. Geometric phase has been 
studied in its entirety encompassing various aspects of evolution including 
adiabatic [7], cyclic, unitary and subsequently the generalization to non- 
adiabatic [8], non-cyclic |9l[l0], and non- unitary [10. .11^ evolution using 



*email:sns@iitk.ac.in 
temail:subhashisli@iitj.ac.in 



1 



either a dynamic or a kinematic or approach |12l 111) . Non-unitary evo- 
lution in the context of open systems has also been reported [13] • These 
studies have been carried out both for pure and mixed states [141115) . On 
the other hand an operational definition of the geometric phase has been 
given by Sjoqvist et al [IB]. There have been further studies on the rela- 
tionship of the geometric phase of an entangled system and its sub-system 
[17) . Recently, a proposal for studying the Pancharatnam phase and its 
relation to non-local quantum correlations has been suggested [TS] in the 
two-photon interferometric set up. 

Holonomic quantum computation and entanglement require the evolu- 
tion of qubits in parameter space by the control of parameters which are 
physically feasible. Both geometric phase and entanglement are impor- 
tant tools in quantum information processing. Atom photon interaction 
provides common ground for exploring the relationship between geomet- 
ric phase and entanglement with the added advantage of the existence of 
control parameters for manipulating the two photon states. In this let- 
ter we study the non-adiabatic, non-cyclic, non-unitary evolution of the 
geometric phase of the two-photon state corresponding to the two modes 
emitted by a three level cascade system interacting with two driving fields 
[19| . The two photon state is in general a mixed state, however it may also 
be prepared in a pure state by proper control. We consider the evolution 
in the state space by varying the control parameters namely, the driving 
field strength and detuning. 

2 The Model 

The scheme considered here corresponds to a three level cascade system 
interacting with two coherent fields which address the only two allowed 
dipole transitions |i)-f->-|i-|-l),i = l,2 with energy separation given by Ui. 
Two counter propagating ( Doppler free geometry) driving fields of nearly 
equal frequencies wli and iJL2 and respective strengths Qi and Q2 are 
resonant with these two transitions. The decay constants of the energy 
levels 1 3) and |2) are indicated by Fa and r2 respectively. The parameters 
Al, A2 refer to the detunings of the driving fields. This scheme may be 
realized, for example, in Rb vapor with the corresponding energy levels 
5Si/2,5P[i/2 and 5-D5/2 and has been used by Banacloche et al. [22] for 
studying electro-magnetically induced transparency. The parameters used 
for the numerical computation correspond to this scheme which we label 
as scheme I. 

The driven three level atomic system emits radiation corresponding to 
the dipole transitions 1) — >■ |i).The resulting field is a two photon state 
in the Fock space and its density matrix is of rank three which follows from 
Schmidt decomposition. This class of two photon states is in a smaller 
sub-space and requires just (A'^-l- 1)^ — 1 variables rather than 4^ variables 
to completely determine the N — 2 state. This reduction simplifies the 
analysis enormously. 

The available control parameters for tailoring the two photon state are 
the driving field strengths Qi, i = 1,2 and the detunings Ai, i = 1, 2. 
corresponding to the two transitions. Let us denote the two photon de- 
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Figure 1: (Color online) Three level cascade system corresponding to the ^'^Rb 
atoms driven by two fields cji and 102- 

tuning by A = Ai + A2. At two photon resonance A = 0. Some of the 
non-classical properties of this two photon field have been reported in the 
literature [201 119) . In particular, a detailed numerical study of the entan- 
glement distribution as a function of the these control parameters has been 
reported in [19] . In order to establish the dependence of entanglement and 
purity on these control parameters we make a small change in the scheme 
used above, namely, we replace the metastable state (Fa = IMHz) by an 
infinitely long hved state i.e, = which we call scheme II. In this case, 
the two photon decoherence given by (F3 -I- Fi)/2 is zero which results in 
a two photon system with no decoherence. This simplifies the analysis 
and (i) the dependence of the entanglement on the control parameters 
and (ii) the relation between geometric phase and entanglement are more 
transparent. The effect of decoherence is included in the real system. 

The two photon density matrix for scheme I (Fig.l) may be determined 
using the quantum tomographic method described in [1^. We briefly 
recapitulate the procedure here. The atom+field pure state, at a time t, 
is written as 



where i denotes the atomic index while n denotes the photon mode index. 
To be more specific, in the binary notation n — n2 ni, where Ui — 0, 1. 
The two photon state may now be obtained by taking the partial trace 
over the atomic indices. It has been further shown [T5], that the two 
photon density matrix at a time t is equivalent to the atomic density 
matrix at a retarded time t — r/c [21] . Utillizing this equivalence we write 
the two-photon normalized density matrix in the basis {|0), |1), |3)} = 



= n); I = 1, 3; n = 0, 1, 2, 3; 



(1) 



{|00),|01),|11)} as 
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Figure 2: (Color online) Variation of the eigenvalues with A. The solid lines cor- 
respond to scheme II while the dashed lines correspond to scheme I. Parameter 
values are a) and b) fii = 6, 5^2 = 6 and c) and d) fii — 3,ft2 — 6. 



= P^it-r/c) (2) 

Where Af is the normalization factor. The solution for the complex coef- 
ficient, <x"{t) may be obtained by solving the Liouville equation., 

mp^ = [H,p^]~'^{r,p^} (3) 

where F is the relaxation matrix and H is the Hamiltonian given by 

H = ^{uiiai + uj2cr2) + hQi{e~^"^^a^ + h.c.) + 

hn2{e-"^^*at + h.c.) (4) 

in the interaction picture and in the rotating wave approximation. The 
Rabi frequency is given by Q,i = —j-fla+i ■ Ei] i = 1,2, corresponding to 
the two driving field strengths. cr+ = ji -I- 1) err = |i) (i + 1|; i = 1,2 
and o-f = |i + 1) (i + 1| — \i) {i\ are the atomic transition operators. 



3 Geometric phase 

The Pancharatnam relative phase for a pure state is defined as 

a(t) = ^rg({*(0)l*W)). (5) 

The two photon state under consideration, on the other hand is in gen- 
eral a mixed state. We follow here the kinematic approach of Tong et 
al. [11] for determining the geometric phase for non-cyclic, non-unitary, 
non-adiabatic evolution of mixed states. We briefly recapitulate their 
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Figure 3: (Color online) Variation of concurrence in the parameter space (A, X). 
(a) Ideal system (scheme II) and (b) Real system (scheme I) with ^2 = 6. 
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procedure here: A general mixed state density matrix p(t) is subject to 
purification by the introduction of an ancilla as 



I'i'it)) ^2^V\{k)\Mi:))'»Wk);te [0,r]. (6) 

k 

The Pancharatnam relative phase given in Eq. ((5]) reduces to the geo- 
metric phase when the parallel transport condition is satisfied. Tong et 
al. [TT] impose the parallel transport condition for the mixed state to be 
{4>k{t)\d/ dt\(j>k{t)) — 0,k = 1...N corresponding to the N eigenstates. In 
the scheme under consideration = 1,2,3. The geometric phase for the 
mixed state, p"{t), satisfying the parallel transport conditions assumes 
the form 

7s(r) = Arg[Y^ ^ Xk{r)\k{0){MO)\<f>k{r) X 

k 

where Xk (t) are the eigenvalues and (j)k (t) are the corresponding eigenvec- 
tors. Here the evolution is considered over the parameter space consisting 
of {A, X} defined in the next section. 

4 Relation between Geometric phase and 
Entanglement 

In order to get a physical understanding of the dependence of the two 
photon state on the free parameters, we consider first, the density matrix 
of the ideal system (scheme II) near two photon resonance which has the 
simple form 

/ ACS^ -SC(l_-ii72iA) \ 
p ^ i ... -AC^S (8) 

V / 

where S = Sin(X) and C = Cos(X). We have defined the dimension- 
less parameters as X = Tan~^(r2i/r22), A — A/y^Q,f + and 721 = 

^2/2^y^^l + ^^^. The density matrix is first order in the parameter A. 
The reason for choosing the state around this point will become clear 
shortly. Consider the spectral resolution of the density matrix 

p = ^A.lV'0{V'»i; i = l,3 (9) 

i 

where are the eigenvalues and {ipi) are the corresponding eigenstates 
of the density matrix. The steady state eigenvalues at A = is given 
by {0,0, 1}, which means it is a pure state. While in the real system 
(scheme I) the largest eigenvalue is always < 1. However, for any given 
ni,n2, purity is still maximum for A = with the largest eigenvalue 
being maximum at this point. A comparison of the eigenvalues for the 
ideal and real system is illustrated in Fig. 2. The role of the metastable 
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state thus, is to introduce a admixture, due to dccoherencc, to an oth- 
erwise pure state. Since the atomic decay constants arc fixed in a real 
system, the only remaining free parameter for controlling purity is the 
two photon detuning A. The state around the point (A = 0,X) may 
be approximated to a pure state when there is no decoherence. This 
is the reason for choosing to study the evolution near two photon reso- 
nance. In the ideal case ( scheme II) the pure state at A = is given 
by = -Sin(X)|00) + Cos(X)|ll); Sin(X) = Ql/^/^T^■ The 

concurrence for this pure state is C = 2Sin(X)Cos(X) and is maximum 
for X = 7r/4 or fii = which we call the bell state regime since this 
corresponds to the Bell state (— 100) + |ll))/\/2. On the other hand con- 
currence is very small when Qi « fl2. The state is separable in this 
region and we call this the separable state regime. A non-perturbative 
evaluation of the distribution of concurrence in the parameter space for 
the system without decoherence is illustrated in Fig. 3a. It is clear that C 
is maximum in the region where the ratio of the field strengths is approx- 
imately one while the state becomes less entangled as the ratio between 
Qi and O2 differs from one. Thus the parameter X which is a function 
of Qi,il2 is the control parameter for entanglement. Having established 
that A and X are the control parameters for purity and entanglement we 
now show that the geometric phase is related to entanglement. 

We evaluate the geometric phase in the space spanned by the pa- 
rameters {A,X} in the neighborhood of xo ~ {0,X) where \i{0,X) = 
0, VX, i = 1,2 and A3(0,X) — 1, VX. In other words we need to eval- 
uate the geometric phase of the pure state a.t x = iS,X + dX) which is 
essentially the Pancharatnam's phase given by 

7«(X) = Arg[(Mxo)\Mx))] (10) 

here ips corresponds to the pure eigenstate and X = Xo+dx= {S,X + dX) 
is any point chosen infinitesimally close to xo- Taylor expanding 7g near 
Xo and retaining up to second order in 5 and dX we obtain 



r2 

r2 • , r<^„/v^2s. 



dX 

7s(x) = Arg[l- — +p5'-i-y2iCos{XfS- 

Z72iCos(X)Sin(X) S dX 
2 ^ 



2 

(11) 



where 



13 = (-(l/8)Cos(X)^(4-H67li-(5 + 87li)Cos(2X)-h 
Cos(4X)) + C^((l + 87li)Cos(2X) + Cos(4X))/16 

and C is the concurrence at (0,X). We have thus, shown that the geo- 
metric phase 79 (x) at x near xo depends on the concurrence C at xo for 
the ideal system (scheme II). 

We next evaluate the derivative f'g of the geometric phase with respect 
to the two photon detuning A which is shown in Fig.4. To be more 
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Figure 4: (Color online) Variation of d 'jg/d A in the neighborhood of (a) A = 
and X ~ 0, (separable state) and (b) A = and X ^ tt/4 (Bell state). 

specific, we evaluate 7^ in the neighborhood of (A = 0,X). We have 
already mentioned that when X = 0, concurrence C = and we have a 
separable state. Let us examine the behavior of geometric phase around 
this point. The rate of change of the geometric phase in the neighborhood 
of (A = 0, X = 0) varies rapidly for small variations (Fig. 4 (a)) 5 along 
the y axis and seems to be unchanged for small variations dX along the x 
axis. On the other hand when X = n/A concurrence C = 1, which is a Bell 
state. Here, the rate of variation of the geometric phase is relatively very 
slow for small perturbation in both S and dX. This is illustrated in Fig. 4b. 
While in the case of C = 1, the total change of the slope, j'g, is from -1.4 
to -2.7 (Fig.4b), in the case of C = it changes from -3.0 to -14.0 (Fig.4a) 
for the same range of values of 5 and dX. Thus the geometric phase seems 
to change very slowly in the vicinity of maximum entanglement. Faster 
the sweep of the geometric phase or larger the rate of change of geometric 
phase implies weaker the entanglement. In other words, stability of the 
geometric phase in any region of the parameter space seems to indicate 
the state has maximum entanglement in this region. In the general mixed 
case, however, it is difficult to show analytically the relation between 
entanglement and geometric phase. We therefore illustrate the relation 
through numerical simulations. 

5 Numerical Results: 

We present below the numerical results for the real (scheme I) system 
shown in Fig.l. 

We evaluate the three eigenvalues and eigenvectors numerically up to 
all orders in the two photon detuning A and determine the geometric 
phase as a weighted sum of the individual phase corresponding to each 
eigenstate. Since we are addressing hyperfine transitions the eigenvalues 
are always non-degenerate in the parameter regimes of interest, in other 
words the hyperfine separations are larger compared to the Rabi frequency 
of the applied fields. 

All the system parameters are rendered dimensionless by scaling with 
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Figure 5: Variation of the geometric phase and its derivative with the detuning 
parameter Ai. Parameter values of a) and b) are fli = ^2 = 6, A2 = c) 
and d) Qi =3,n2 = 6, A2 = 0, e) and f) r^i = 6, 1^2 = 3, A2 = 0. g) and h) 
ni=n2 = 6, A2 = 3, i) and j) Qi = 1.5, 1^2 6, A2 = 



-0.5 -0.25 0.25 0.5 

Figure 6: Stability of 7g with respect to fluctuations in A and fli — ^l2- 

the atomic lifetime 7 « IMHz. We restrict our study to the parameter 
regime A,Q, < 6.0 which is the lifetime of the first excited state. We 
would like to mention here that in the numerical program, 7^ has been 
initialized to zero for minimum values of the parameter Ai. 

In the previous section, for the scheme I, we have seen that the state 
is dominantly entangled when Qi « Q2 and separable otherwise. Even 
in the case of the real system (scheme I), we may identify three regions 
in the parameter space: (i) when Qi « f22,Ai,2 = the two photon 
state is dominantly in a pure state with very weak entanglement , (ii) 
when r^i ~ f22, Ai,2 = the state is almost pure Bell state of the type 
[|00) — |ll)]/\/2 and (iii) when » f22,Ai,2 = the state is mixed 
with weak entanglement. Non-zero detuning (A 7^ 0) further affects both 
purity and entanglement. A distribution of concurrence as a function 
of the two photon detuning A and Qi — 02 is shown in Fig.Sb for the 
real system with decoherence. The state seems to be show maximum 
entanglement near f^i ~ 0,2 and shows weaker entanglement for ili « O2 
and fii >> ^2. It is clear that even in the presence of decoherence the 
qualitative features of entanglement distribution is similar even though 
there are quantitative changes. This is not very surprising given the fact 
that there is no unique entanglement measure for mixed states. 

We illustrate in detail the variation of the geometric phase in the 
neighborhood of Ai = with A2 = 0. This is effectively varying the two 
photon detuning A. Figs. 5. a, b show the variation of 7g and the rate 
of change of 7^, respectively with Ai near Ai = for fli = Q2 = 6. The 
main feature we observe is that the geometric phase varies slowly near 
the region of entanglement which we showed to be true in the case of the 
ideal system earlier (Fig.Sb). Observe that 7^ is constant in the region 
—0.25 < Ai < 0.25 (Fig.5 a,b). This is substantiated by the vanishing 
of the derivative in this region thus indicating that the geometric phase 
in this case is stable under small perturbations in the neighborhood of 
Ai = 0. The numerical values reveal [T5] that the population is almost 
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equal in the states |00) and with non-zero correlation which is typical 
of a Bell like state. Again in Figs. 5 g, h which corresponds to the 
parameter values fli = Q2 = G and A2 = 3, the geometric phase 7^ 
does not change much and the derivative is nearly zero in the region 
—0.5 < Ai < 0. The entanglement distribution in these two cases is nearly 
the same with the two photon state being an almost pure Bell state. We 
have verified that the variation of the density matrix elements in these two 
cases is less than 2%. In the case of Figs. 5 (c), (d), the parameter values 
are Qi = 3, n2 = 6, A2 = 0. Observe that rate of variation of 7g is more 
rapid compared to (b) and (h) indicating an overall weaker entanglement. 
The slope decreases to a minimum at Ai — 0, implying that the state is 
maximally entangled here. However, 7g is never constant in this regime 
which means it is not stable under perturbation in the parameter A. A 
look at the density matrix elements reveals [ID] that the state is an almost 
pure state but with weak entanglement. 

Lastly, we address the question whether geometric phase can distin- 
guish between pure separable state and mixed state. As shown if Fig. 3b, 
the distribution of concurrence for the real system does not seem to dis- 
tinguish between the pure separable state and the mixed state. In fact 
this is true of other measures like negativity and fidelity. However, the 
variation of geometric phase seems to show this distinction. For the pa- 
rameter values fli — 6,Q2 = 3 the state is a mixed state [19] and the 
variation is illustrated in Fig. 5 e,f. Here 7g is never a constant and the 
slope which corresponds to the rate of change of 7^ shows sharp fluctu- 
ations in the neighborhood of Ai — 0. While for a pure separable state 
r^i — 1.5, Q2 = 6, the derivative shows a smooth variation as illustrated 
in Fig. 5 i, j. Thus geometric phase is unique in the sense that it carries 
information on both purity /mixedness and entanglement. 

In Fig. 6 we present the details of this plateau where 7g is stable. The 
plateau in the center marked zero corresponds to the region where 7g does 
not change. This is sorrounded by region where 7g changes by 1%. The 
negative sign indicates the decrease in 7g. The figure illustrates a total 
decrease of 7g by 10%. Thus, in a given interval of the control parameter, 
faster the sweep of 7g (with respect to the control parameter), weaker the 
entanglement and vice versa. This is consistent with the distribution of 
concurrence shown in Fig. 3b. 

6 Summary and Conclusions 

We have studied the evolution of the geometric phase in the general set 
up namely non-adiabatic, non-cyclic and non-unitary evolution of mixed 
states. The scheme corresponds to a Hamiltonian system consisting of in- 
teracting atom-photon system. We have identified the control parameters 
for manipulating the purity /mixedness and entanglement of the emitted 
two-photon states. We have been able to show the relation between the 
geometric phase and entanglement in terms of the control parameters an- 
alytically for an ideal system without decoherence. For the real system, 
we show that the variation of the geometric phase is correlated to the vari- 
ation of the concurrence on the control parameters Jli and Ai, i — 1,2. 
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We have shown that 7g is robust under fluctuations of these control pa- 
rameters only for pure Bell like states. We further show that the variation 
of the geometric phase is able to distinguish between pure separable and 
mixed states. These features provide a diagnostic tool for determining 
the entanglement of the two-photon state. It is interesting to note that 
such a connection, between entanglement and geometric phase, has also 
been found, recently [23], in another system, viz. Heisenberg spin chains 
undergoing unitary evolution. 
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